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GaAs-based two-dimensional electron gases (2DEGs) show a wealth of remarkable electronic 
states 1 3 , and serve as the basis for fast transistors, research on electrons in nanostructures 4 5 , 
and prototypes of quantum-computing schemes 6 . All these uses depend on the extremely 
low levels of disorder in GaAs 2DEGs, with low-temperature mean free paths ranging from 
microns to hundreds of microns 7 . Here we study how disorder affects the spatial struc- 
ture of electron transport by imaging electron flow in three different GaAs/AlGaAs 2DEGs, 
whose mobilities range over an order of magnitude. As expected, electrons flow along nar- 
row branches that we find remain straight over a distance roughly proportional to the mean 
free path. We also observe two unanticipated phenomena in high-mobility samples. In our 
highest-mobility sample we observe an almost complete absence of sharp impurity or defect 
scattering, indicated by the complete suppression of quantum coherent interference fringes. 
Also, branched flow through the chaotic potential of a high-mobility sample remains stable 
to significant changes to the initial conditions of injected electrons. 

Scanning gate microscopy (SGM) images of electron flow in two-dimensional electron gases 
(2DEGs) 8 18 provide direct spatial information not available in conventional electrical transport 
measurements. Our SGM studies show how varying disorder affects electron flow, and enable us to 
infer information about the disorder potential in our different samples. Achieving a detailed picture 
of the disorder potentiaii2l22! ma y help to understand why exotic electron organization emerges in 
some 2DEGs and not others, and to aim for ever weaker disorder or even tailored disorder^. 
By analyzing the differences between images of flow in our samples, we find that the higher- 
mobility samples are increasingly dominated by small-angle scattering instead of hard- scattering^. 
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Finally, we investigate an unusual property of electron flow through such a small-angle scattering 
disorder potential: though the disorder potential is classically chaotic, branches of flow are stable 
to significant changes in initial conditions. 

On each of three 2DEG samples defined in GaAs/AlGaAs heterostructures (Table [T]), we use 
a home-built scanning gate microscope to image the flow of electrons emanating from a split-gate 
quantum point contact (QPCj^HIH at 4.2 K, as schematically shown in Fig. Id. Using a recently 
established techniquePES we measure the conductance across the QPC while scanning a sharp 
conducting tip ~20 nm above the surface of the sample. We negatively bias the tip to create a 
depletion region in the 2DEG below. When the tip is above a region of high electron flow from the 
QPC, it backscatters electrons through the QPC, reducing the measured conductance. By scanning 
the tip and recording the drop in conductance AG for each tip location, we thus image electron 
flow. Images of flow gathered in this way have been found to accurately reproduce the underlying, 
unperturbed flow patterns^^l Unless otherwise stated, the conductance of the QPC is set to the 
first plateau, 2e 2 /h. Samples with mobilities as high as samples B and C or as low as sample A 
have not been previously imaged using this technique. 

Even in relatively clean 2DEGs, in which the mean free path is much greater than the Fermi 
wavelength, electrons do not propagate smoothly. Instead they often flow along caustics - narrow 
branches of concentrated current density - which form owing to the cumulative effect of passing 
over many small-angle scattering sites that focus electrons into these branches^ 10 ! 25 ! 26 !. The disorder 
potential is composed of a low, bumpy background of small-angle scattering sites from ionization 
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of some of the Si donors present above any conventional GaAs-based 2DEG, as well as sharp, 
tall spikes, causing hard-scattering, from unintentional impurity atoms embedded in the 2DEG, 
more prevalent in lower- mobility 2DEGs 22 . The length scale of potential fluctuations in the 2DEG 
produced by ionized donors in the donor layer is set by the separation between the 2DEG and the 
donor layer, a distance which is significantly shorter than the mean free path. 

Figure 1 shows the striking differences in flow patterns for electrons traveling in 2DEGs 
with a wide range of mobility and mean free path. The flow varies from twisted and diffusive in 
the lowest-mobility sample (sample A) to straight, smooth branches of flow in the highest-mobility 
sample (sample C). Here we measure that the characteristic length l B over which branches remain 
straight increases in samples with longer mean free path I, as might be expected. As a simple 
way to estimate Ib, we average the distances along branches between all observable pairs of points 
where one branch splits into two, in five images each for samples B and C (see Supplementary 
Information, Section III); we identify 49 and 29 branch points in all images for samples B and 
C, respectively. We find in sample B with I = 13 //m that Is = 480 nm, and in sample C with 
I = 28 /im that l B = 740 nm. 

We now use SGM to estimate the density of hard-scatterers in our highest-mobility samples, 
something which is difficult to do from basic transport measurements. In Fig. 2, we discover an 
almost complete absence of hard- scattering from impurities and defects in the highest-mobility ma- 
terial by comparing interference fringes which decorate images of flow in the three samples. All 
previous sufficiently high-resolution images of electron flow through QPCs using this backscat- 
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tering imaging technique have shown interference fringe J^EIl. However, as shown in Fig. 2, in 
our highest-mobility sample, sample C, the interference pattern is completely suppressed. Sample 
A, with the lowest mobility, shows strong interference fringes throughout the flow pattern, while 
sample B, with intermediate mobility, shows interference fringes only in a few locations. 

As schematically shown in insets to Fig. 2, and as detailed in previous worliML-y 25 ! 27 !, our 
imaging technique yields fringes due to the interference of two sets of paths- path 1 (blue), the 
roundtrip path electrons take from the QPC to the tip and back, and paths 2 (red), the sum of all 
roundtrip paths electrons take from the QPC to all fixed scatterers in the sample and back. When 
we move the tip away from the QPC by half the Fermi wavelength, Xp/2, the phases of paths 2 
do not change, but the phase accumulated along path 1 increases by 2ir; thus, we expect a fringe 
spacing of Xp/2. As has been shown previously 10 , the distance from the QPC over which these 
fringes persist is not cut off by thermal broadening (see Supplementary Information, Section VI) 
and can be expected to be as long as the quantum coherence length, which is 9 //m for sample B 
and 6 fim for sample C at 4.2 KPl 

In high-mobility sample C (Fig. 2c), we observe no interference fringes because the disorder 
potential contains very few sites that hard backscatter electrons back through the QPC: paths 2 
apparently have much smaller amplitude than path 1. In sample B, interference fringes are only 
visible in a few rare locations, one of which is shown in Fig. 2b, signifying the presence of a 
hard-scatterer near the center of interference fringes. The rare appearance of fringes is consistent 
with this high-mobility sample containing a sparse sprinkling of impurities close to the 2DEG. 
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Judging from the density of observed scattering centers in the flow through five devices in sample 
B (see Supplementary Information, Section IV), we roughly calculate the density of hard-scatterers 
is ~ 3//im 2 , which is close to estimates in similar materials^. In sample C, the lack of any 
observable fringes indicates that the density of hard-scatterers is less than ~ 0.5//im 2 . Sample A 
contains too many scattering sites to establish a reliable estimate. 

Figure 3 shows an interesting consequence of small-angle scattering in a high-mobility sam- 
ple (sample B): the branched flow pattern is robust to changes in initial conditions. Previous ex- 
periments have inferred information about electron trajectories in quantum chaotic systems from 
transport measurement^, but here we use SGM to directly image these trajectories and their sen- 
sitivity to initial conditions. Building on previous experiment^ 12 ! 13 !, we change the initial position 
and angle at which electrons are launched into the 2DEG by applying unequal potentials to the two 
QPC gates. By putting as much as a 500 mV offset between the two gate voltages, simulations 
indicate (see Methods section) that we can shift the position of injected electrons by ±30 nm, and 
the initial injection angle by ±5°, away from the more negatively biased gate. A 60 nm shift is 
significant in this system because it is similar to both the 50 nm correlation length of our simulated 
disorder potential and the ~65 nm width we calculate for the QPC channel. 

Electron trajectories in a low, bumpy disorder potential are classically chaotic: very large 
changes in the flow pattern might be expected from this significant change in initial conditions. 
Contrary to these expectations, the experimental flow pattern appears to be remarkably immune 
to changes in electron injection position and angle. As we shift the QPC potential, electrons flow 
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along the same set of branches, with the only difference being changes in the relative strength of 
flow along each branch. Figure 3a-c shows experimental images for which the gates are unevenly 
charged and electrons from the first conductance plateau are pushed to the left, not pushed and 
pushed to the right, respectively. Comparing the images in Fig. 3a-c, it is remarkable that the 
central branches 4 fxm away from the center of the QPC (visible in Fig. 3a-d, denoted with brown 
arrows) deviate by less than 10 nm, despite our estimate that the center of the QPC moves 60 nm. 
Figure 3d shows electron flow when the QPC is symmetrically opened to the second conductance 
plateau, 4e 2 /h, so that a wider spread of electrons emanates (through both the first and second 
conductance modes). The positions of the left branches (visible in Fig. 3a,d, denoted with a blue 
arrow) and right branches (visible in Fig. 3c,d, denoted with a green arrow) are extremely stable 
as well. 

To understand these data, in Fig. 4 we numerically simulate flow in a low, bumpy disor- 
der potential chosen to match the physical characteristics (Fermi energy 7.5 meV, mean free path 
13 jum, and donor setback 25 nm) of sample B (see the Methods section). In Fig. 4a,b, we show the 
chaotic nature of the potential by displaying the extreme sensitivity of individual classical trajec- 
tories (shown in red) to small changes in initial conditions. We calculate that the Lyapunov length, 
the average length scale over which two trajectories initially separated by a small displacement 
exponentially diverge, is 300 nm for this potential (see Supplementary Information, Section VIII). 
Based on these results, we might expect branch shapes and positions on the scale of these im- 
ages to be highly sensitive to initial conditions. Contrary to this expectation, the overall simulated 
branched flow pattern (shown in blue) in Fig. 4a and Fig. 4b, composed of an ensemble of many 
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classical trajectories with an angular spread similar to that expected from a QPC, demonstrates 
remarkable immunity to changes in initial conditions. This stability is in accord with previous 
theoretical work on stability in chaotic systems which found that while individual trajectories 
are unstable the overall flow pattern, or manifold, can be much more stable to variations in initial 
conditions. 

However, whereas classical mechanics has previously been found adequate to explain branch 
formation 10 and here reproduces stability of branches to small changes in initial position, it cannot 
reproduce the full experimentally observed stability of flow patterns. The simulated classical flow 
pattern in Fig. 4c, which is shifted by 60 nm, is very different from that in Fig. 4a, demonstrating 
that the branches in classical flow are not stable to a full 60 nm shift in initial position of a point 
source (see Supplementary Information, Section II for more information). We find that a quantum 
treatment is necessary to reproduce the observed experimental robustness, as in Fig. 4d,e. The 
positions of most of the branches in the simulated quantum flow pattern deviate by less than 10 nm 
even several microns away from the QPC, even though there is less than 5% overlap between the 
injected electron \ip\ 2 in the two simulations (see inset Fig. 4e). We posit that not only diffraction 
out of the QPC but also quantum scattering off the disorder potential steadily increases the spatial 
overlap between the two differently injected beams as they move away from the QPC, providing 
an additional source of stability in the quantum simulations. 

In conclusion, the technique we describe here, which enables our probe to act as a detector 
for individual hard- scattering sites, could prove useful as a measurement of the impurity concentra- 
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tion in the highest- mobility 2DEGs, and could provide important insight into the limiting factor of 
the ultimate mobility in GaAs-based 2DEGs. Also, the branch stability we observe could be impor- 
tant in applications which require guiding electrons through 2DEGs and related low-dimensional 
semiconductor structures. 

Methods 

To determine the potential produced at the 2DEG by gates on the surface of the sample, we use 
SETE, a three-dimensional Poisson-Schrodinger solver written by M. Stopa (Director of National 
Nanotechnology Infrastructure Network at Harvard University). SETE first determines the pro- 
file of the conduction band edge along the growth direction of the hetero structure by solving the 
Schrodinger equation in one dimension and then calculates the potential and electron density at 
the 2DEG self-consistently from the Poisson equation^. We then quantum-mechanically time- 
evolve the first mode of conductance through the QPC potential computed by SETE. Based on the 
resulting quantum-mechanical wave function we estimate how much we shift the initial position 
and angle of injected electrons by unevenly charging the QPC gates. In addition, we determine an 
appropriate spread of angles to use with our classical trajectories: the full- width at half-maximum 
of electron trajectories emanates at ±10°. 

To simulate electron flow in sample B, we create the disorder potential from a simple phys- 
ical model based on two experimentally known parameters: the distance between the 2DEG and 
donor layer, and the mean free path (see Supplementary Information, Section VII, for more infor- 
mation). We consider only donors in the donor layer, and assume each ionized donor is metallically 
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screened by the 2DEG, creating a density perturbation corresponding to a potential dip with a full- 
width at half-maximum of 40 nm. We convolve a random donor distribution with this potential. 
This results in a disorder potential with correlation length 50 nm. We then set the amplitude of the 
potential energy fluctuations so that the quantum-mechanically simulated mean free path in this 
disorder potential matches the measured value in our real sample. We use this disorder potential in 
our classical simulations, and in the quantum-mechanical simulations we add it to the gate poten- 
tial determined by SETE. In the quantum-mechanical simulations, we then quantum-mechanically 
time-evolve a wavepacket with energy spread equal to 3.5 kT to account for the effects of thermal 
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Figure 1 Electron flow in samples with different mobilities, a-c, Electron flow originates from the 
QPC formed between the two depletion regions, schematically indicated in black below each image. The 
length over which branches remain straight depends on the sample's mean free path, a, Sample A, with 
mean free path 1.5 fim, has branches that quickly change direction or are hard to assign a direction, b, 
Sample B, with mean free path 13 fim, has branches that remain straight over moderate length scales. 
The average distance along a branch of flow between observable branch points is 480 nm. c, Sample C, 
with mean free path 28 /im, has branches that remain straight over the longest length scales. The average 
distance along a branch of flow between branch points is 740 nm. b,c, Red arrows denote two neighboring, 
observable branch points where one branch splits into two. d, Schematic showing metallic gates and tip 
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(orange) creating depletion regions (black) in the 2DEG (green) buried below the surface of the sample. 
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Figure 2 Interference fringes in samples with different mobilities, a-c, Different samples show vary- 
ing strength of interference fringes, correlated with the samples' mobilities. The presence of interference 
fringes more than the thermal length away from the QPC indicates the presence of hard backscattering 
in the sample due to impurities or defects, a, Interference fringes in sample A representative of fringes 
decorating all images of flow in this sample, b, Electron flow in sample B, which shows interference fringes 
in some portions and not in others. Visible interference fringes are rare in sample B; this panel focuses on 
one of these rare regions. Fringe spacing is denoted by red circles to the left of electron flow, c, Electron 
flow in sample C representative of the absence of interference fringes in this sample. Insets: Two classes of 
paths interfere to form fringes (see text): path 1 , blue; paths 2, red. Black regions indicate depletion areas 
from the gates and the black disc indicates depletion from the tip. Gray areas indicate hard-scattering sites 
in the sample, due to impurities close to the 2DEG for example. The lack of hard-scattering sites (gray dots) 
in the higher-mobility samples (B and C) suppresses interference fringes in these images. 
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First plateau I Second plateau 




Figure 3 SGM images of branch stability in sample B. This figure shows the remarkable experimental 
stability of branched flow. The reference grid (200 nm) denotes the same physical location for precise 
comparison of electron flow between images, a-c, Images of electron flow with the QPC open to the 
first plateau of conductance and electrostatically shifted to the left, not shifted and shifted to the right, 
respectively, d, The second plateau of conductance with gates symmetrically energized. Black regions 
indicate schematically the depletion regions forming the QPC. The thick red arrows denote the initial average 
position and momentum of injected electrons. On changing the initial conditions of injected electrons, the 
location of branches does not change, but the current density flowing along each branch does. In all images, 
the central branches, one of which is labeled with a brown arrow, are visible. In a and d the left branches, 
labeled with a blue arrow, are visible. In c and d the right branches, labeled with a green arrow, are visible. 

19 



I Classical simulations 1 Quantum simulations 




Figure 4 Classical and quantum-mechanical numerical investigations of branch stability. These 
figures show simulations of electron flow into a disorder potential chosen to resemble that expected in 
sample B. See Supplementary Information, Section I, for log-scale plots. The reference grid (200 nm) 
denotes the same physical location in every simulation for precise comparison of electron flow between 
simulations, a-c, Classical electron trajectories are injected from a point source with 19° angular spread. 
The point source in b is 2 nm to the right of that in a, and the point source in c is 60 nm to the right of 
that in a. The trajectory initially pointing directly up in a and b is highlighted in red. The slight shift in initial 
conditions between a and b causes this trajectory (red) to deviate significantly after only 1-2 ^m of travel, 
but the overall branched flow (blue) formed in a and b still looks surprisingly similar. However, with a full 
shift of 60 nm (the calculated shift between Fig. 3a and c) the branches formed in c do not resemble those 
in a. d,e, Quantum-mechanical simulations of flow through an electrostatically shifted QPC into the same 
disorder potential as a-c. The injection points in d and e correspond to those in a and c, respectively. Most 
branch patterns in d and e remain stable between the two simulations. Some similar branches are denoted 
with dashed dark green lines. Inset in e: Lateral profiles of electron probability density (IV'I 2 ) immediately 
after injection in d and e. Electrons are injected 60 nm further to the right in e than in d owing to uneven 
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voltage on the gates (for comparison, X F = 56 nm). The similarity in flow patterns is significant and striking 
because the overlap between the two injection profiles is minimal. 
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Sample 


A 


B 


C 


density, n (10 11 cm -2 ) 


4.3 


2.1 


1.5 


mobility, \i (10 6 cm 2 /Vs) 


0.14 


1.7 


4.4 


mean free path, I (/im) 


1.5 


13 


28 


distance from surface to 2DEG (nm) 


57 


68 


100 


distance from donors to 2DEG (nm) 


22 


25 


68 



Table 1 : Properties of 2DEG samples, measured in a Hall bar configuration. Sample A was grown by the 
commercial grower IQE. Sample B was grown by H. Shtrikman. Sample C was grown by L. N. Pfeiffer and 
K. W. West. Each of the three heterostructures hosts a 2DEG at a single GaAs/AIGaAs heterojunction, not 
in a square quantum well. 
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